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Abstract 

Floquet theory is used to describe the unstable spectrum at large scales of the /?- 
plane equation linearized about Rossby waves. Base flows consisting of one to three 
Rossby wave are considered analytically using continued fractions and the method 
of multiple scales, while base flow with more than three Rossby waves are studied 
numerically. It is demonstrated that the mechanism for instability changes from 
inflectional to triad resonance at an 0(1) transition Rhines number Rh, independent 
of the Reynolds number. For a single Rossby wave base flow, the critical Reynolds 
number Re"^ for instability is found in various limits. In the limits Rh — > oo and 
A; ^ 0, the classical value Re'^ = \pl is recovered. For Rh and all orientations 
of the Rossby wave except zonal and meridional, the base flow is unstable for all 
Reynolds numbers; a zonal Rossby wave is stable, while a meridional Rossby wave 
has critical Reynolds number Re'^ = For more isotropic base flows consisting 
of many Rossby waves (up to forty), the most unstable mode is purely zonal for 
2 < Rh < oo and is nearly zonal for Rh = 1/2, where the transition Rhines number 
is again 0(1), independent of the Reynolds number and consistent with a change 
in the mechanism for instability from inflectional to triad resonance. 
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1 Introduction and the main results 



The /3-plane equation is a model for large-scale geophysical flows, in which 
the effect of the earth's sphericity is modeled by a linear variation /3 of the 
Coriolis parameter. Due to the j3 term, solutions of the linearized /3-plane 
equation are dispersive waves, called Rossby waves, with anisotropic dispersion 
relation given by <^(k) = —pk^/k"^, where <^(k) denotes the wave frequency. 
Zero-frequency linear eigenmodes are zonal flows and form a slow manifold. 
In the absence of the P term, the model reduces to the two-dimensional (2D), 
incompressible Navier Stokes equations. 

The present work is motivated in part by numerical simulations of /3-plane 
flow driven by stochastic forcing at small scales (Chekhlov, Orszag, Galperin 
and Staroselsky [1] , Smith and Waleffe [2] , Marcus, Kundu and Lee [3] , Huang, 
Galperin and Sukoriansky [4] ) . These authors show that the presence of the f3 
term leads to anisotropic transfer of energy to large scales and accumulation 
of energy on the slow manifold corresponding to zonal flows. The generation 
of slow, large scales is also characteristic of other models of geophysical flows, 
including three-dimensional (3D) rotating flow [2], and the 2D and 3D Boussi- 
nesq equations for rotating stably stratified flows [5], [6]. The accumulation 
of energy in slow modes is interesting given that resonant triad interactions 
cannot transfer energy from fast waves directly to slow modes with zero fre- 
quency as shown, for example, in Longuet-Higgins and Gill [7] for the /9-plane 
equation. Newell [8] showed that higher-order interactions, such as resonant 
quartets, can lead to energy transfer from Rossby waves to zonal flows. To- 
wards a complete understanding of the simulations [1], [2], [3], [4], this paper 
describes the large-scale instability of base flows maintained by a single or 
several deterministic force (s), for values of P ranging from < (3 < oo, and 
for flnite Reynolds numbers. 

The results herein extend previous stability analyses of the 2D Navier Stokes 
equations and the /3-plane equation. Using Floquet theory and a continued 
fraction formulation, Meshalkin and Sinai [9] studied the linear stability of the 
2D Navier Stokes equations with a Kolmogorov base flow, and found the crit- 
ical Reynolds number Re'^ = \f2 for large-scale instability with wavenumbers 
k For Reynolds numbers slightly above the critical value, the instability 
occurs for modes with wavevector k perpendicular to the wavevector of the 
Kolmogorov base flow. Friedlander and Howard [10], Belenkaya, Friedlander 
and Yudovich [11] and Li [12] used similar techniques to study the stability of 
the 2D inviscid Navier Stokes equations (the 2D Euler equations). Sivashin- 
sky [13] used Multiple Scales analysis on the perturbation equation for the 2D 
Navier Stokes equations with a Kolmogorov base flow, and obtained a non- 
linear amplitude equation describing large-scale flow with Reynolds number 
slightly above the critical value Re'^ — -\/2. Generalizing the work of Sivashin- 
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sky [13] to the /?-planc, nonlinear amplitude equations for super-critical, large- 
scale flow were derived by Frisch et al. [14] for a zonal Kolmogorov base flow 
and small /3 <C 1, and by Manfroi and Young [15] for a nearly meridional 
Rossby wave base flow and flxed P — 0(1). Instability of resonant triad inter- 
actions is absent from [14] because the base flow is zonal and /3 is small, and 
from [15] partly because the Reynolds number is only slightly supercritical 
(slightly larger than -\/2, see 5). 

In all of the above stability analyses, the large-scale instabilities have a pre- 
ferred direction because the base flows are anisotropic (consisting of a single 
mode). More isotropic base flows for the 2D Navier Stokes equations were con- 
sidered by Sivashinsky and Yakhot [16], who found that large-scale instability 
is not observed for any finite Reynolds number if the base fiow is sufficiently 
isotropic. They used Multiple Scales analysis for the linearized perturbation 
equation, with isotropic scaling in space and periodic base flows consisting 
of several modes. Recently, Novikov and Papanicolau [17] have extended the 
linear analysis of [16] to nonlinear analysis for a rectangular base flow. 

For inviscid /?-plane flow with large Rhines number, Lorenz [18] found growth 
rates of linearly unstable modes from a low-order truncation of the inflnite 
dimensional eigenvalue problem associated with a base flow consisting of a 
Rossby wave and a constant zonal flow. Gill [19] considered a Rossby wave 
base flow, and using the continued fraction formulation, noted the transition 
from inflectional instability for Rh — > cxd to resonant triad instability for 
Rh^O. 

In this paper, using the continued fraction formulation and a Rossby wave base 
flow with wavevector p, we obtain the equation for the linear growth rate at 
the response wavevector k. We study the equation at various Rhines numbers 
Rh oc 1//3 to flnd the growth rate analytically or numerically, and obtain 
critical Reynolds numbers in different limits. For Rh — > oo and k 0, the 
critical Reynolds number approaches the value for the 2D Navier Stokes 
equations. Manfroi and Young [15] have also recently studied linear stability 
of the /3-plane equation in the Rh — > oo limit, but with Kolmogorov base 
flows instead of Rossby wave base flows. Using Multiple Scales analysis in the 
limits Rh oo and A; ^ 0, they find a critical Reynolds number smaller than 
the value a/2 for the 2D Navier Stokes equations. Our analysis is compared to 
that of Manfroi and Young [20] in Section 5. For small Rhines numbers (large 
P), we discuss the instability of resonant triad interactions. We also generalize 
the work of Sivashinsky and Yakhot [16] for more isotropic base flows to the 
/3-plane equation. 

The organization of the paper and the main results are as follows. In Sec- 
tion 2, the linear stability problem is formulated. For forcing of a single 
mode, Floquet theory and continued fractions lead to an inflnite dimensional 
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eigenvalue problem. For finite Reynolds numbers, convergence of the contin- 
ued fractions is guaranteed, and the convergence rate is increased for large- 
scale instability with A; <C 1. (Convergence of the continued fractions is dis- 
cussed in Appendix A.l.) In Section 3, we consider the base flow wavevector 
P — iPx,Py) — (cos(q;p), sin(Q;p)) and response wavevector k, with k <^ 1 and 
fixed Rh. In this case the stability problem can be well approximated by a 
low dimensional (3 x 3) eigenvalue problem that may be solved analytically. 
This truncated system is expanded in powers of k to find the growth rate of 
the perturbation with A; <^ 1. For fixed Re and Rh, this growth rate is given 

by 



^{a{k; ap, Re, Rh)) = - 



Re 



1 - 



Re'' 



2 {l + Re^u'^y- 



■X 



(1 



8cos^a + i?eV)] + 0{k^) 



where ^{z) denotes the real part of a complex number z, ci;(k) is the wave 
frequency given by the dispersion relation cj(k) = —Rh'^k^/k'^, and a is the 
angle between the base flow wavevector p and the perturbation wavevector 
k. The marginal Reynolds number for long- wave instability Re'^ (k; ap, Rh) 
corresponds to zero growth rate, and is found by setting to zero the lowest- 
order term in the above expression for the growth rate 3?((7(k; ccp, i?e, 
For fixed ctp and Rh, Re^(k; ap, Rh) for k ^ is the same as the critical 
Reynolds number for A; — > 0. We show that Re^^kiap, Rh) is always greater 
than the critical Reynolds number Re'^ — y/2 for isotropic, 2D, Kolmogorov 
flow with k ± p [9], [16], and [21] (see also [14] and [15]). Indeed, for Rh ^ oo 
we find 



lim i?e°(k; ap, Rh) = ^2, k ± p. 

R/i— >oo 



For fixed Rh and |k| — > 0, we find 



lim mm Re^(k:a^, Rh) 

ko-*0\k\<ko ^ ' 



= < 



cos^ Q;p(l — 8 sin^ ap) 



tana 



Ph 



< Clip < sin ^ J 



sin-i I < ftp < 7r/2 



which is a generalization of the result for the critical Reynolds number Re'^ — 
y/2 for 2D isotropic fiow in the limit k ^ [16]. 



In the next Section 4, we consider |k| < 1 and 0.1 < Rh < oo by numerical 
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computation of growth rates for a less severe truncation (31 by 31) of the 
infinite dimensional eigenvalue problem. An important observation is that the 
nature of the instability to large scales changes from an inflectional instabil- 
ity to a triad resonance as the Rhines number is decreased from Rh — oo 
(Figure 4), where the transition Rhines number is 0(1) (See Gill [19]). For 
forcing wavevector with ap ^ 7r/2 and Rh < 1, there is a band of unstable 
wavevectors near k resonant with ±p, and the width of this band decreases 
as Rh decreases. As shown in reference [7] , resonant triad interactions cannot 
transfer energy to modes with = 0, and indeed the numerical computations 
show the most unstable wavevector k near the resonant trace has 7^ 0. For 
large Rhines numbers and Reynolds numbers above critical {e.g., Re — 10), 
the inflectional instability is strongest for ap — and perturbation wavevector 
with kx = 0. 

The hmit Rh — > is studied in Section 5. The main result is that for Rh — > 0, 
the critical Reynolds number is given by 

V2, ftp = 
0, < ttp < 7r/2 
00, ccp = 7r/2. 

The critical Reynolds number is 00 for ap = 7r/2 in the sense that for any 
Reynolds number, the base flow is stable if Rh is small enough. Thus, resonant 
triad interactions reduce Re'^ below the value y/2 for small Rh and < ctp < 
7r/2. 

In Section 6, we compare our results with the recent work by Manfroi and 
Young [15] using Multiple Scales analysis. In [15], Manfroi and Young consider 

the /9-plane with Kolmogorov base flows, rather than the Rossby wave base 
flows considered herein. Because their choice of base flow is different from 
ours, the critical Reynolds number for /c — > and Rh — > cxd is different from 
\/2. We confirm their results, and reproduce our results for Rossby wave base 
flows using Multiple Scales analysis. 

In Section 7, base flows consisting of several Rossby waves are considered. 
The growth rate of large-scale instability is found analytically for large Rh 
following [21]. In Appendix A. 3, the method of [21] in Fourier space is shown 
to be equivalent to the Multiple Scales analysis following [16]. An important 
result is that an "isotropic" base flow consisting of three Rossby waves is 
unstable above Re > ^32/3, in contrast to the case of pure 2D flow, which is 
stable for a similar base flow [16]. We also numerically studied the growth rate 
of large-scale instability for base flows of many Rossby waves (up to forty), 
as a function of Re and Rh. For 2 < Rh <^ 00 and large enough Re, the 



lim Re^{ap, Rh) = < 

Rh — >0 
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most unstable large-scale mode is piircly zonal with kx = 0. For small Rh < 1 
and large enough Re, the most unstable large-scale mode has wavevector k 
with kx small but nonzero and kx < ky, again indicating a tendency towards 
large-scale zonal flows. A summary is given in Section 8. 



2 Linear stability of the /3-plane equation 

2. 1 Base flows and perturbations 



The governing equation for the stream function ( V x ^'z = u) in the /3-plane 
approximation is given by 



(V^^)* + J(V2^, ^) + i3^x = i^V^^ + (V X F) • z (1) 

where J {A, B) — AxBy — AyBx is the Jacobian, (3 is the linear variation of the 
Coriolis parameter, v is the viscosity and F is a force. For any nonzero vector 
a, a = a/|a| is the unit vector in the direction of a. With the length scale 
and the velocity scale U°, the corresponding nondimensional form of (1) is 



(V^*)* + J(V2*, + = ^V^^ + ^(V X F) • z (2) 

Rn Re Re 

where the Reynolds number Re and the Rhines number Rh are given by 



Re — — and Rh — — -— . 

vp 13 

The inviscid hmit of (2) with F = yields wave solutions, called Rossby waves, 

U(x, t- k) = (k X z)t/(k) exp(«^(x, t; k)) + c.c. 
where c.c is the complex conjugate and the phase ^(x, t] k) is given by 

^(x,i;k) =k-x-a;(k)i (3) 

with the dispersion relation a;(k) 
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In this paper each wavevector k = (kx, ky) has the polar coordinates represen- 
tation (A;, ctk) so that k = A;(cosQ;k,sinQ;k). 

Consider a base flow U°(x, t) for (2) given by a sum of m Rossby waves of 
unit wavevectors Pj — (cos ap. , sin ap. ) , 



U°(x,t;{p,}™0 = ^(p, X z)sin^(x,i;p,). (5) 

i=i 

Throughout this paper, we always assume that wavevectors Pj are unit vectors 
unless otherwise specified. This fiow has a stream function 



^r°(x,i;{p,.}- J ^-^cose(^,t;pj) (6) 

3=1 

and can be maintained by forcings at m wavevectors Pj (and their conjugates) 
with the following force: 



F(x, t; {p,}™ i) = J2iPj X ^) sin^(x, t; p,). 

Note that *° satisfies V^*° + = 0. In particular, for m = 1 and Rh ^ oo 
{P — 0), this fiow is a sinusoidal Kolmogorov fiow for the 2D Navier Stokes 
equations. 

A small perturbation t) = ^(x, t) — \I'°(x, t) satisfies the following lin- 
earized perturbation equation. 



(VV)t + + V'V)y - ^^(V^ + V'^), + - i?e- W = 0, (7) 

where the quadratic nonlinear terms in ip have been dropped for the linear 
stability analysis. To study the linear stability at wavevector k, we use Floquet 
theory to seek a solution ip of the form 



V'(x, t; {p,}- J = eV^-'''"^^ Yl V'n exp ^ n,-^(x, t; p,) j , (8) 

where the first sum is over n = {rii, • • • , n^). Inserting this form of ip into (7), 
the following infinite dimensional linear eigenvalue problem is obtained 



A -I- it!e q (n) + iojn ■0: 
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-j^ rn 

I^(Pj X q(n - %))(g^(n - Gj) - 1)V'„- 



m 

where q(n) = k + Z^^i^jPj, g(n) = |q(n)| and e^- is the j-th unit vector 
(Cj = (0, • • • , 1, • • • , 0)). Also, is given by 



m 

= -uj{k) - I 5]nja;(pj) | +uj{q{n)). (10) 



Our goal is to find the eigenvalue A of (9) with the largest real part for each k 
with fixed Re, Rh and pj, and this eigenvalue is denoted by o"(k; {ap^}^i, Re, Rh). 
Thus, the linear growth rate is given by its real part 3?((T(k; {ap^.}^^. Re, Rh)) 
where 'R{z) is the real part of any complex number z. If the growth rate is 
positive for k, the base flow is unstable to a perturbation with wavevector k. 

Although it is not easy to analyze (9) for general m, the case of m = 1 can 
be studied using the technique of continued fractions following Meshalkin and 
Sinai [9]. For m — 1 with p = pi, the base flow (6) is 

^'°(x,i;p) = -cos^(x,i;p), (11) 



and (9) reduces to 



[A + i?e q (n) + iix)n]il^n - 



p X k • z 

2g2(n) 



(g^(n-l)-l)^„_i 



P''''-^g2(n+l)-l)V'n+i = 0. 



2q^{n) 



[12) 



Here, q(n) = k + np, q{n) — |q(n)| and 



Un = —^i)^) - nuj{Y)) + a;(k + p). 



(13) 



We define an(A, k, a^. Re, Rh) and (in(k, a-p, ipn) by 
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(14) 



Then, (12) can be written as 

(J-ndn + dji-l — dn+i = 0. 

Formally, one can show that a„ satisfies 



11 

oo + T + 1 = 0. (15) 

ai + „ . ' 1 a-i + „ . ' 1 

"■2+ „„.... °-2+ „ 



This idea was first explored by Meshalkin and Sinai [9] for the 2D Navier 
Stokes equations in the special case when k and p are perpendicular to each 
other. They obtained the classical critical Reynolds number for the first 
instability at large scales. Later, Friedlander et al. [10], Belenkaya et al. [11] 
and Li [12] used similar techniques to study the stabihty of a Kolmogorov flow 
for the inviscid (i?e = oo) Navier Stokes equations (Eulcr equations). In this 
paper, we consider only the viscous [Re < oo) case to have fast convergence of 
the continued fractions in equation (15). The convergence is straightforward 
since for fixed A, k, ctk, -Re and we have 

|a„| ^ cx) as |n| ^ oo 

due to the viscous term Rer^q{nY in the numerator of a„. In Appendix A.l, 
we discuss the convergence of the continued fractions in equation (15). For 
each A satisfying (15), one can construct {ipn} so that (8) is well-defined. For 
this discussion, we refer to [10], [11] and [12]. Unfortunately, the technique of 
continued fractions can not be extended for m > 1. 



2.2 Numerical computation of linear growth rates 



For all figures of numerically computed linear growth rates, contours of the real 
part of (T(k; {ap^ j^^. Re, Rh) are plotted. To indicate the unstable wavevec- 
tors, the shaded regions consist of k with 3?((T(k; {ctp^}^!. Re, Rh)) > and 
the darker regions consist of k with 

3?((7(k;{Q;p.}^i,i?e,i?/i)) > 0.9 max ^{a{\i]{ap.}f^^,Re,Rh)) 
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where the maximum is taken over all k under consideration in each figure. 
Thus, the wavevector with the strongest instability is inside the darker regions. 

For base flow (11) with m — 1, the linear growth rate 3?(a(k; ctp, i?e, 
is the real part of the solution to (15) with the largest real part. Since we 
consider the viscous case {Re < oo), the continued fractions in (15) have fast 
convergence. Thus, we can truncate the continued fractions in (15) for \n\ > N 
for some large A^. This is equivalent to truncating the infinite dimensional 
eigenvalue problem (12) for \n\ > N and considering the {2N + 1) x (2N + 1) 
finite dimensional eigenvalue problem. We find all eigenvalues for the truncated 
eigenvalue problem using the MATLAB function 'eig'. Then, (j(k; Op, Re, Rh) 
is the eigenvalue with the largest real part. For most numerical results, we 
use N = 15. We check the results with = 16 for some cases and find no 
significant difference between the results for N ~ 15 and N — 16. 

For the real part of (T(k; ccp. Re, Rh), one can check the symmetry 

3fJ((T(k; Op, Re, Rh)) = ^{a{-k; ctp. Re, Rh)) 

from the form of our perturbation (8) with m = 1. This enables us to con- 
sider wavevectors in the half-plane. Another interesting observation is that 
K((T(k; ttp. Re, Rh)) is invariant under the shift by jp for j = 0, ±1, ±2, • • •. 
By the form of our perturbation (8), one can also show that if A is an eigenvalue 
of (12) for k = ko, then 

A + i(-a;(ko) - ju{p) + u(ko + jp)) 

is an eigenvalue of (12) for k = ko-|-jp. Thus, 3?((T(k; ap. Re, Rh)) — ^{a{'k.+ 
jp; Op, Re, Rh)) for any j = 0, ±1, ±2, • • •. As an example to illustrate these 
properties. Figure 1 shows contours of ^{a(k; ap. Re, Rh)) for ap = 7r/4, 
Re = 10 and Rh = 2 in -2 < < 2 and -2 < ky < 2 with a grid 
of size 200 x 200. From Figure 1, we observe the symmetry and the shift 
invariance. Thus, it is sufficient to consider wavevectors k with |k-p| < 1/2 and 
(p X k) ■ z > 0. For special cases ap = and Op = 7r/2, we have the additional 
symmetry ^{a{{kx, ky); ap, Re, Rh)) = ^{a{{±kx, ^ky); ap. Re, Rh)) which is 
due to the symmetry of a;(k) = —Rh'^k^/k^- This allows us to consider only 
wavectors k with < k^ < 1/2 and ky > for ctp = 0, and A;^ > and 
< ky < 1/2 for CKp = 7r/2. It is not easy to find the growth rate numerically 
for m > 1 since the construction of the matrix for the eigenvalue problem of 
(9) is nontrivial. Thus, later in Section 7, we use "severe" truncation of (9) to 
find approximate growth rates for m > 1. 
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Fig. 1. Contour plot of K(cj(k; Op, Re, Rh)) for for CKp = 7r/4 at Re = 10 and Rh = 2. 
3?((T(k; ctp, i2e, Rh)) is positive in the shaded regions. 



3 Large Scale Analysis for fixed Rh 



In this section, we derive the growth rate for small A; <^ 1 with base flow 
(11) at fixed ap, Re and Rh by asymptotic expansion of equation (15) of 
continued fractions. From the growth rate, we also find the critical Reynolds 
number in the limit |k| ^ 0. Prom (15), we show that the linear growth rate 
3?((7(k; CKp, Re, Rh)) is given by 



3f?((7(k; ap, Re, Rh)) = - — 

Re 



Re' 



sin^ a 



Re'iu 



2\2 



(l - 8 cos^ a + Re^uj"^) 



(16) 



where u = CL^(k) is the dispersion relation (4) and a = — is the angle 
between the response wavevector k and the base fiow wavevector p. The higher 
order terms H3 and H4 satisfy 



for some nonnegative constants C3 and C4 depending on ap, Re and Rh but 
not on k. The exphcit expression for is given below (see (29)). In particular, 
C3 is proportional to Rh~^ and is zero when Rh = 00. Thus, if Rh = 00, then 
the growth rate (16) is reduced to 



K(cr(k; ctp, Re, Rh = oo)) 




sin^Q;(l -Scos^a) + 0{k^) 



which is the large-scale growth rate obtained by Sivashinsky and Yakhot [16] 
and by Waleffe [21]. Note that no assumption is made on Rh for (16). However, 
since the higher order term H-^ is 0{k^ / Rh), the leading order term in (16) 
is a good approximation of the actual growth rate if Rh is large. The leading 
order term can also be derived using Multiple Scales analysis by assuming 
that the Rhines number is large as shown in Appendix A. 3. See [6], [16], [14], 
[15] and [20] for more details. Note that one advantage of using equation (15) 
of continued fractions is that we can obtain the growth rate (16) up to third 
order in k explicitly. 

3. 1 Derivation of the growth rate 

This section is devoted to deriving the growth rate formula (16) from (15) for 
A; <^ 1. We will assume that 



This assumption is required as follows. In general, the infinite dimensional 
eigenvalue problem (12) has infinitely many eigenvalues for each wavevector 
k. In particular for k = 0, (12) has eigenvalues 

0, -1^1 Re, -2^ /Re, 

and thus 'U{a{Q;ap,Re,Rh)) ~ 0. Therefore, for /c <g 1 we can expect that 
(12) has eigenvalues 

0{k), -1^1 Re + 0{k), -2Vi?e + 0{k), 

Assumption (17) guarantees that we obtain the first eigenvalue which has 
the largest real part for /c ^ 1. In fact, with assumption (17) one finds that 
3f?(cr(k; ttp, i?e, i?/i)) = 0{k'^) as shown below (see (27)). 

For A; ^ 1, we show that (15) can be reduced to 



3?(cT(k;ap,i?e,i?/i)) > 



1 



(17) 



2Re 




(18) 



The remainder R^ satisfies 



12 



|i?5| < C^k', 

where C5 is a constant depending on ap and Re but not on Rh and k, and 



lim — - = 0. 

Re^o Re 



To this end, we have 



a{k]ap,Re,Rh) + {1/Re)q^{n) +iujn > ^ n = ±1, ±2, . . . (19) 

oRe 

by assumption (17). The expression inside the absolute value is a part of the 
numerator of a„ in (14). Other expressions in (14) satisfy 



pxk- z = 0{k), q^{±l)-l = 0{k), g^(±n) - 1 = 0(1), n = 2,3,..., 
for /c <^ 1. Using these and (19), one can show that a„ in (14) satisfy 



^ <CRek^ = 0{Rek^) 



(20) 

^ <CRek = 0{Rek), for n = ±2, ±3, . . . 



for some constant C independent of Rh. Using (20), one finds that 



^ ^0{Rek), (21) 



so that we have 



— + 0{Re^k^). (22) 



The details of the derivation of (21) and (22) are given in Appendix A.l. 
Equation (18) follows from (15) using (21) and (22). 

The growth rate (16) is obtained from (18) by asymptotic expansions of a±i 
in k, as will now be shown. A Taylor expansion of (13) gives 
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u!±i — — c<;(k) ^ Ci;(p) ± a;(p ± k) 



: -a;(k) + Du;{p) ■ k ± -li^D^u;(p) ■ k + 0{k^/Rh) 



(23) 



where Dlo{^ and D^cj(p) are the first and second derivatives of u; at p re- 
spectively and k^ is the transpose of k. From (23) one can see that uj^ and uj^i 
are the same up to 0{k'^/Rh) for A; ^ 1. For simphcity, let Vti = Duj{p) ■ k, 

^2 = ^k^D^uj{p) ■ k and a = a{k; ap, Re, Rh). Note that = 0{k/Rh) and 

^2 = 0{k^/Rh). From (14), oq and a±i can be rewritten as 



a±i 



2k{a + k^Re) 
sin a{l — k'^) 



2(1 ± 2A; cos a + P) (o- + (1 ± 2k cos a + P) /i?e + ia;±i) 



(24) 



sin q;(±2 cos a + k) 



For simphcity, we set 



q± = Vl±2k 



cos a 



do = (J + Re 



-1 



(i±i = ±2A; cos a/Re + ifli , (i±2 = k'^/Re^ ifl2 
Expanding l/a±i in A; and using (23) and 25, we find that 



(25) 



1 A;^ sin a (±2 cos a + A;) 

~ 2ql^ {do + cii + 4 + 0{k^/Rh)} 

k"^ sin a(2 cos a ± k) 

2ql,do(l + ^ + ^ + 0{kyRh) 
\ do do 



k"^ sin a 
2 do 



(±2 cos q; + A;) (l ^ 2A; cos a + A;^(4 cos^ a - 1) + 



k"^ sin a 



2do 

2iQi cos a 



±2 cos a + /c < 1 - 4 cos^ a 1 + 



do 



do 



±Gi + 1G2 



+ 0{k'), 



(26) 



where Gi and 6*2 are given by 
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^ 2 COS a ( ^ , ^ 2kcosaVtil Re\ kQi 
Gi = — ; — {Q2 + 2kcosafli-\ ' ' 



do 



1 



do J do ' 



G2 = ^ I 2 cos a(il2 + 2A; cos — A;f2i + , 
do \ do 



Assumption (17) guarantees that the denominator do = a + Re ^ — icu is not 
zero. 

Using (26) and (25), (18) leads to 



a — 



sin^ a 



+i 



Re 2{a + Re'^ - iu) 
k sin^ a 



1-4 cos a 1 + 



Re-^ 



a + Re~^ — iu 



2(a + Re-^ - iu 



-G2 + 0(r). 



(27) 



Prom expression (27), we see that |(7(k; ctp. Re, Rh) \ — 0{k'^/Re) as mentioned 
before. Expanding the denominator {a + Re~^ — icu), we obtain 



<7(k) = - 



Re 



1 - 



sin^ a 



2(i?e-i - iu) 



1-4 cos a 1 + 



Re 



-1 



Re"^ — iuj 



fcsin^o; ~ ^/,4x 

+'7^r^i — + 0(k^) 
2[Re — icu) 



(28) 



where G2 is given by 



1 



Re ^ — iu 



2 cos Q;(Jl2 + 2k cos aili) — k^x + 



8A; cos^ Q;i7i/i?e' 
i?e~^ — iu 



Note that is Oik^jRh) since Qi = 0{k/Rh) and ^2 = 0{k'^/Rh). In 
particular, (72 = if it!/i = 00. Finally, the growth rate (16) is obtained by 
taking the real part of a in (28). The third order term in (16) is given by 



\ 2[Re ^ — lu) . 



(29) 



and H4 is the fourth order remainder. Note that is 0(k^/Rh) since G2 — 
0{k^/Rh). In particular ^ ii Rh ^ 00. 
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3.2 Marginal Reynolds numbers i?e°(k; Op, Rh) 



The critical Reynolds number is the smallest nonnegative Reynolds above 
which a base flow base becomes unstable. With base flow (11), the critical 
Reynolds number depends on ap and Rh, and is denoted by Re'^{ap, Rh). 
Similarly, the critical Reynolds number i?e^(k; ctp, Rh) at a fixed wavevec- 
tor k is defined by the smallest nonnegative Reynolds number above which 
^{a(k] ap, Re, Rh)) is positive, and can be found by setting the growth rate 
equal to zero. In other words, if the Reynolds number is larger than i?e'^(k; Op, Rh), 
then the wavevector k becomes unstable. If such a Reynolds number does not 
exist at k, we define Re'^(k; ap, Rh) = oo in the sense that the wavevector k 
is stable for all finite Reynolds numbers. Note that Re'^(k; ap, Rh) = oo does 
not imply that the wavevector k is unstable for the inviscid {Re = oo) /3-plane 
equation. In terms of Re'^(k; ap, Rh), Re'^^ap] Rh) satisfies 

Re'^{ap, Rh) = niini?e'^(k; ap,Rh). 



The critical Reynolds number i?e'^(k; ap, Rh) for A; <^ 1 may be found by 
setting the growth rate (16) to be zero. However, it is nontrivial to find 
Re'^(k; ap, Rh) from the growth rate (16) since the explicit expression for the 
higher-order term H^j^ in (16) is not known. Instead, we study the smallest non- 
negative Reynolds number which make the leading order term of (16) equal 
to zero by solving 



„ o 2(l + Re'^ujY , 9 9 9 . ^ 

Re^ = — 2 — — ^ — -, and 1 - 8 cos^ a + Re^u^ > 0. (30) 

sin a |1 — 8 cos^ a + ite^'o;^} 

It is called the marginal Reynolds number and denoted by Re'^ik; ap, Rh). 
Similar to the critical Reynolds number, we define i?e°(k; ap, Rh) = cxd if no 
Reynolds number satisfying (30) exists for k. In particular, for Rh = oo, the 
marginal Reynolds number is given by 



i?e°(k; ap, Rh = oo) = J — 5 — — - for 1 — 8 cos^ a > 0, 

y sin a(l — 8 cos^ a) 

and its minimum is equal to \/2 for a = 7r/2 as shown by Sivashinsky and 
Yakhot [16] and by Waleffe [21]. In general, it is clear that Re^(k; ap, Rh) 
is different from the critical Reynolds number Re'^(k;ap, Rh) but can be un- 
derstood as the critical Reynolds number in the limit — > O.Q Since the 



Re^(k;ap, Rh) is the critical Reynolds number in the limit /c — > in the sense 
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higher order term in the growth rate (16) is proportional to 1/Rh, Re^ is 
meaningful for large Rhines numbers. 



We show that the marginal Reynolds number satisfies 



i?e°(k;ap,i?/i) > ^2, 



(31) 



where the value \/2 is the classical critical Reynolds number for the 2D Navier 
Stokes equations [9], [16], and [21] (see also [14] and [15]). This inequality can 
be easily shown as follows: if (1 — 8cos^ a + Re'^u'^) > 0, then from (30) 



i?e°(k; ap,Rh) 



Re'^uj^ + 2(2 - sin^ a)Re^uj'^ + 2(1 - sin^ a(l - Scos^ a)) 
s\v? q;((1 - 8 cos^ a) + Re^uj'^) 



> 



since (2 — sin^ a) and (^1 — sin^ a{l — 8 cos^ a)^ are nonnegative for all a. In 
particular, when a — — ap = t:/2, we have 



lim 

Rh—*oo 



i?e°(k; CKp, Rh) 



2. 



(32) 



Then, (31) and (32) imply that for any ojp. 



lim min Re°{k; ap, Rh) = V2, (33) 

Rh^oo k 

recovering the classical critical Reynolds number \/2 for the 2D Navier Stokes 
equation. This is different from the recent result by Manfroi and Young [20] 
for the /3-plane equation with a steady Kolmogorov base flow. They observe 
that the critical Reynolds number is less than y/2 in the limit Rh oo. We 
compare their results with ours in Section 6. 



3.3 Critical Reynolds numbers in the limit k ^ 



In this section, we consider the marginal Reynolds number i?e°(k; ctp, Rh) for 
A; — > 0. Without loss of generahty, we only consider < ctp < 7r/2 and the 



that 



lim min Re^(\s.:ar,,Rh) = lim min Re'^(\s.:ar,,Rh). 
feo-»0|k|<feo fco-*0|k|<jfco 

This can be shown by an argument similar to the one in Section 3.3. 
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result is the following: 



lim min i?e*^(k: 

feo^O|k|<fco 



cos^ ap{l — 8 sin ctp) 



, < Op < sin 



-1 1 



(34) 



8| tanttpl, 



sin~^ I < CKp < 7r/2. 



Note that (34) depends only on Op but not Rh and is an increasing function in 
< ctp < 7r/2. For the case ctp = tt/2, (34) is infinite in the sense that for any 
Reynolds number, wavevectors k with |k| < ko are stable if ko is small enough. 
By definition, the critical Reynolds number Re'^{ap, Rh) is less than or equal to 
(34). The expression (34) increases to infinity as ctp increases to 7r/2. Thus, we 
expect that Re'^{ap, Rh) is strictly less than (34) if ctp is not too close to zero. 
In this case, the first unstable wavevector is observed "away" from the origin 
as the Reynolds number is increased from zero. For example, for ctp = 7r/4 
and Rh — 10, instability is first observed near wavevector k = (—0.17,0.18) 
at Reynolds number Re ^ 1.59 (Figure 2). We have 



Re%a = 7r/4, Rh = 10) < 1.59 < 8 = lim min i?e"(k; ap = 7r/4, Rh = IC 

feo-^O |k|<A;o 

As the Reynolds number is increased, the unstable region becomes a long 
strip and reaches to the origin at the Reynolds number (34) (Figure 2(b) with 
ctp = 7r/4 and Re = 8). This is different from the 2D Navier-Stokes equations 
with a Kolmogorov base flow where the flrst instability is always observed at 
/c — at Reynolds number ^/2 [9], [16], and [21]. 

The rest of this section is devoted to deriving (34). Since ^{a{k; ap, Re, Rh)) — 
3?(a(--k; ap. Re, Rh)) as discussed in Section 2.2, we only consider < < tt. 
For a fixed real number r], we introduce a function 



fi'^'^v) = 7;Tr, — {l - 8cos^77 + 'u;| for w > 0. 
2 (1 + w)^ ^ ' 

Then, f{Re'^uj'^; a) is the reciprocal of the right-hand side of (30). Recall that 
uj — u!{k.) — —Rh~^kx/k^ and consider all wavevectors k = A;(cosQ;k,sinQ;k) 
with A; > and a fixed a^- In the case Rh < oo and ap ^ 7r/2, the range of 
a;(k) consists of all positive numbers, and in the case «k = 7r/2, the range of 
a;(k) is zero. In any case, the minimum of i?e°(k; ap, Rh) over wavevectors k 
with a fixed cck satisfies 
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(a) Re = 1.67 



0.25 



(b) Re = 8tan7r/4 = 8 

X 10"' 



0.15 




0.05 



Fig. 2. Contour plot of ^{a{]^■,ap,Re,Rh)) for ap = 7r/4 at Rh = 10 and (a) 
Re = 1.67 (b) Re = 8. In (b), note that the x-axis is scaled by 10~^ and the y-axis 
is scaled by 10~^. 



mm 

fc>0 



Re^{k{cosap, sin ap); ap, Rh) 
The maximum of f{w; rj) is 



> 



1 



max^>o f{w]a)' 



(35) 



max /(w; 77) = < 

tu>0 



and achieved at 



I sin^ 77(1 — 8 cos^ 77), for cos^ ^7 — 
^ tan^ 77, for cos^ 



(36) 



i(;(77) = < 



0, 



for cos^ V Tf 
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16cos^77 — 1, for cos^77 > 



(37) 



To show the equahty in (34), we demonstrate that the inequality "<" and ">" 
both hold. The inequality "<" is trivial for ap — 7r/2 since (34) is infinite. If 
ap < sin~^(l/4), then we have 



min (i?e°(k; ap, Rh)f < min (i?e°((0, ky); ap, Rh)f 



|k|<feo 



\ky\<ko 



cos^ Clip (1 — sin^ dip) 



since c<;(k) =Q\ikx — Q- Next, we suppose that ap > sin (1/4). Motivated by 
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(37), fix the Reynolds number Re — Stanccp and consider the curve consisting 
of k's satisfying 



Re'^uj(k)'^ — 16 cosset — 1, kx < 0, ky > 0, a = cck — ccp. 



(38) 



Since we have a;(k) = —Rh'^kx/k'^-, the curve (38) is tangential to the hne 
/c-r = at the origin (actually, this curve is along the unstable region in 
Figure 2(b).). At a fixed wavevector k on curve (38) with Re — 8tanQ;p, we 
have 



Re 



Re^ sin^ a 
2 (1 + Re^uj^Y 



{l - 8cos^Q; + i?eV} 



Re 



2 tan^ a 
1 — 64 tan a, 



P 64 



Re 



tan (ftk — cip) 
tan^(7r/2 — ctp) 



> 



since we have 7r/2 < at < tt on the curve. This imphes that i?e°(k; ctp, Rh) < 
8tanci;p and establishes the inequality "<" in (34). 

Now, consider the opposite inequality ">" in (34). For this purpose, let {k" = 
A;"(cosQ;kn, sincKkn)} be a sequence satisfying 



(a) lim |k"| = 0, 



< (b) lim min Re^(\s.:aj,, Rh) = lim i?e°(k", -R/i), 

fco^O|k|<fco ^ ' ^ n-»oo V ' P' /' 



(c) Re^(k^;ap, Rh) is increasing. 
For example, {k"} can be chosen by the relation 



(39) 



min Re%k; ap, Rh) = Re%k''; ap, Rh) and |k"| < n 

|k|<n~i 

We suppose that 



-1 



lim ctk" 



TT 



(40) 



and show this at the end of this section. From (40) we have 



lim — 7r/2 — a^, 



(41) 
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where — ctk" — ctp. Applying (39a), (35), (36) and (41) in order, we have 



hm min i?e°(k; ctp, Rhf = hm (i?e°(k", Op, Rh)f 

fco— >0 |k|<feo n— >oo 



> hm 



n^oo max^>o/(w;Q;") 
2 



hm < 

n— »oo 



sin^ q;"(1 - Scos^ a") 
64 



, if cos a" < J 



if cos a" > I 



= < 



— 5 -, if sinojp < J 

cos^ Q;p(l — 8 sin ctp) 



if sinojp > i 



64tan^ ctp, 

which estabhshes the inequahty ">" in (34). 

We are left to show (40). Suppose (40) does not hold so that we have 



lim cosctk" ^ 0. 

n— »oo 

Then, since /c'^ — >• by (39a), we have 



lim |a;(k")| = lim 



1 cos Ok" 



oo. 



Rh k" 

By definition, M„ = i?e°(k"; ctp, RhY satisfies 



(42) 



sin^ a" 



2 (1 + M„a;(k")2)2 
Using (31), (39c), (42) and (43), we have 



{l -8cosa'* + M„cu(k")2} = 0. 



(43) 



= lim 

n— »oo 



1 - 



Mn sin^ a" 



2 (1 + M„u;(k")2)2 
This is a contradiction and thus (40) holds 



{l -8cosQ;" + M„a;(k")2} 



1. 
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4 Linear Stability as Rh decreases 

Prom the numerical solution of (15), we consider the behavior of the growth 
rate of the perturbation at wavevectors k over |k| < 1 for various Rhines num- 
bers Rh. We first give a preliminary discussion of resonant triad interactions 
between wave solutions of the /3-plane equation. 

4-1 Resonant Triad Interactions 

As stated in Section 2, the /3-plane equation has wave solutions, called Rossby 
waves, 

*(x, t; k) = ^(t; k) exp (i6l(x, t; k)) + c.c, 

where c.c is the complex conjugate and the phase ^^(x, t; k) is given by (3). We 
may represent the solution ^'(x, t) in terms of the waves by 

*(x,i) = ^^(i;k)exp(i^(x,i;k)) with ^*(i;k) = ^{t;-k) 

k 

where * is the complex conjugate. Then, the governing equation (2) with F = 
is written as 

(a. + i^.^)^(,k) = 

J2 Ckpq^*{t; p)^*{t; q) exp (z(cu(k) + cu(p) + u;{q))t) (44) 

k+p+q=0 

where the coefficients are given by Ckpq = {q^ — P^)[(p x q) ■ z]/(2A;^). In this 
section p is not necessarily a unit vector. These coefficients satisfy 

k'^Ckpq + P^Cpqk + q^Cqkp — 0, 

k^Ckpq + p'^Cpqk + q^Cqkp ~ 

as can be deduced directly from energy and enstrophy conservation by triad 
interactions. See Waleffe [22], [23] for more discussion. From (44), in the limit 
Rh — > 0, energy transfer between three wavevectors k, p and q is maximal 
when the frequency ci;(k) + ci;(p) + ci;(q) = 0. Otherwise, a nonzero frequency 
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^■-1 ^05 05 1 

Fig. 3. Resonant traces with ap = 0,7r/4 and 7r/2 

produces an oscillation and energy transfer between the wavevectors will be 
zero on average. 

A triad (k, p, q) is resonant if it satisfies the resonant conditions 

k + p + q = 0, ij{k) +cu{p) +uj{q) = 0. (45) 

The dispersion relation (4) of the /5-plane equation allows resonant triads as 
shown in Longuet-Higgins and Gill [7]. They showed that for a given wavevec- 
tor p with px 7^ 0, the trace of wavevectors resonant to p forms a closed 
curve. For p with = 0, the resonant trace consists of the wavevectors k 
with ky — —py/2. Figure 3 shows three traces resonant to p for ccp = 0, 7r/4 
and 7r/2 with p = 1. For < ap < it/2, as illustrated in Figure 3, resonant 
traces are tangential to the line k^. = at the origin. By asymptotic expansion 
of (45) in kx and ky for k ^ 1, one finds that 



(46) 

which implies that 



k^^O{¥), ak = 7r/2 + 0(A;'-^) (47) 

where r = 3 for < ctp < 7r/2 and r = 4 for ftp = 0. The resonant trace for p 
with a-p = 7r/2 is different from ^ ti/2 and it consists of the modes with 
ky = — |. In this case, if the triad (k, p, q) is a resonant triad, then k — q 
which results in Cpqk = from (44). Thus, the mode p with Px = can not 
receive or lose energy by direct resonant triad interactions ([7], [22]). 
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(a) Re = W, Rh = oo (h) Re = 10, Rh = 2 (c) Re = 10, Rh = 1/2 




(d) Re = 10, Rh= 1/10 (e) Re = WOO, Rh = 2 {{) Re = 1000, Rh = 1/2 




Fig. 4. Contour plot of 5R((T(k; ap, i?e, for Op = 0. Resonant traces including 
ibp are solid lines for (c), (d) and (f). 

4-2 Growth rates for various Rh and |k| < 1 



We study the linear growth rate over k with |k| < 1 at various values of Rhines 
numbers by solving the eigenvalue problem (12) numerically as discussed in 
Section 2.2. Equation (12) is truncated at = 15, resulting in a 31 x 31 matrix. 
The Reynolds number is finite and large enough to exhibit instabilities. In most 
figures, the Reynolds number is fixed at 10 unless otherwise specified. 

First, we consider the case = in Figure 4. Recall that growth rate 
3fJ((j(k; ttp = 0, Re, Rh)) has symmetries about both the a;-axis and the y-axis 
(Section 2.2). Thus, it is enough to consider wavevectors k with < A;^^ < 1 
and < ky < 1 in each plot of Figure 4. We use a grid size of 100 x 100. At 
Rh = oo (Figure 4(a)), the most unstable wavevectors are along the line of 
kx = which is perpendicular to the base fiow wavevector p, consistent with 
the results by Sivashinsky and Yakhot [16] and by Waleffe [21]. We refer to the 
instability at wavevectors near perpendicular to p at large Rhines numbers 
{Rh > 2) as the infiectional instability. This terminology is motivated by the 
existence of infiection points in the profile of the sinusoidal base fiow, which 
can be considered as the source of instability.^ 

As Rh is decreased, the maximum growth rate decreases (Figure 8) and the 
^ By Rayleigh's Inflection Point Criterion. See [24] for example. 
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(a) i?e = 10, = 2 (b) i?e = 10, = 1/2 




Fig. 5. Contour plot of ^{a{k.] Op, Re, Rh)) for ctp = 7r/4 at Re = 10. Resonant 
traces including itp are solid lines in (b). 

most unstable wavevector is not exactly at = but near = (Fig- 
ure 4(c)). As Rh is decreased more, the unstable region becomes much nar- 
rower (Figure 4(d)). In fact, this narrow band lies along the traces resonant 
to the base flow wavevector p and its conjugate —p. The resonant traces cor- 
respond to black lines in Figure 4(c), (d) and (f). Thus, the nature of the 
instability changes from the inflectional instability to the resonant triad inter- 
action as Rh is decreased. This transition of instabihty is observed at Rhines 
number between 2 and 1/2 for Reynolds number Re — 10. Actually, the or- 
der of the Rhines number of the transition appears to be independent of the 
Reynolds number. Figure 4(e) and (f) exhibit the transition at Rhines number 
between 2 and 1/2 with Reynolds number Re = 1000. See also Gill [19]. 

A similar transition is observed for any < < 7r/2 as demonstrated in Fig- 
ure 5 for ftp = 7r/4 at Re = 10 over — 1 < A;^; < 1 and < ky < 1 with a grid 
size of 200 X 100. As Rh is decreased from 2 to 1/2, the most unstable wavevec- 
tor occurs near the traces resonant to ±p. Thus, the most unstable wavevector 
moves from a wavevector perpendicular to p (an inflectional instability) to a 
wavevector resonant to ±p. 

The transition of instability to the resonant triad interaction can be shown 
from (15) for Rh — > as follows. Suppose that the wavevector k is not resonant 
to either ±p so that uj±i 7^ in (14). Then, since |a;±i| — > 00 in the Rh — > 
limit, we have |a±i| — > 00 in (14). Thus, (15) is reduced to 

oo = 0, 
which is equivalent to 

a{k; ttp. Re, Rh) = -k'^/Re < 0. 

Thus only wavevectors resonant to the base flow wavevector p (or its conju- 
gate) can be unstable. 
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(a) Re = 10, Rh^2 



(b) i?e = 10, Rh^ 1/2 




Fig. 6. Contour plot of 3fi((T(k; ap,Re, i?)) for ap = 7r/2 at Re = 10. 



The case of = 7r/2 must be considered separately. As discussed in Section 
4.1, the base flow wavevector p = (0, 1) cannot lose energy directly by resonant 
triad interactions, and as the Rhines number is decreased, the flow becomes 
stable as show in Figure 6 at Rh — 1/2. Indeed, at any finite Reynolds number, 
the zonal base flow {ap = vr/2) becomes stable at sufficiently small Rhines 
number. Stability of zonal flows for small Rhines number was shown by Kuo 
[25] and is discussed further in Section 5. 

The transition from inflectional to resonant triad instability can be observed 
with only two triads:k, p, — (k + p) and k, — p, — (k — p). This corresponds to 
the severe truncation of (15) to 



1 1 

ao + — + = 

Oi a_i 



(48) 



which is same as the eigenvalue problem obtained from (12) after 3x3 trunca- 
tion with n — —1, 0, 1. Figure 7 shows the hnear growth rate 5R((j(k; ctp. Re, Rh)) 
by solving (48) and the change in the nature of the instability is observed 
as Rh is decreases. In particular, the Rhines number for the transition for 
Q!p = is also between 2 and 1/2, as in Figure 4. Because of the severe 
truncation, 3?(cr(k; ctp. Re, R)) is not invariant under the shift jp. Comparing 
Figure 4(b), (c) and Figure 7, at least qualitatively, (48) is a good approxima- 
tion of (15) for wavevectors k with k < 1. This observation is the background 
of the numerical computation in Section 7 to solve (9) approximately with 
m > 2. Note that such truncation (3 x 3) was also used in Lorenz [18] and 
Gill [19]. 
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(a) Re - 10, Rh^2 (b) Re = 10, Rh = 1/2 




Fig. 7. Contour plot of 5R(o"(k; Op, i?e, obtained from (48) for ap = at 
Re = 10. Resonant traces including ztp are solid lines in (b). 

4-3 Inflectional instability with k ± p 

As defined in section 4.2, we refer to an inflectional instability as the instabil- 
ity at wavevectors near perpendicular to p at large Rhines number {Rh > 2). 
Inflectional instability with k ± p is observed to be strong for large Rh. For 
example, for ctp = and Rh > 2, inflectional instabilities are dominant (Fig- 
ure 4). For «p = 71 /A and Rh > 2, the wavevector with the strongest growth 
rate is almost perpendicular to p, although they are not exactly perpendicular 
to p (Figure 5). Thus, studying the inflectional instability with k ± p can give 
a qualitative picture of the behavior of instability for large Rhines number. 

Figure 9 shows contours of growth rates 9?((T(k; ap, Re, Rh)) due to the inflec- 
tional instability with k _L p at Reynolds number Re = 10. For each k with 
< kx < 1 and < ky < 1, p is chosen so that ccp = cck — 7r/2 (k _L p). For 
Rh — oo, the growth rate is invariant under rotations since the Navier Stokes 
equations are isotropic (Figure 9(a)). For Rh = 2, note that the maximum 
growth rate at Re = 10 is strongest for ctp = and monotonically decreasing 
for < ctp < 7r/2 (Figure 8). Thus, the inflectional instabilities with k _L p for 
Rh — 2 are strong for k with k^ — {ap — 0) in Figure 9(b) since the most 
unstable wavevectors are near perpendicular to p for large Rh. For Rh — 1/2, 
wavevectors k resonant to ±p have strong growth rates and such wavevec- 
tors can not have kx = (Section 4.1 and 4.2). In Figure 9(c) the strongest 
inflectional instability with k _L p is observed at wavevectors k with kx ^ 0. 

We conclude that for large Rhines numbers {Rh > 2 for example), the inflec- 
tional instability is the dominant mechanism for instability (Section 4.2) and 
the strongest inflectional instability is observed when ap = and ak = tt/2. 
This observation is consistent with the picture of anisotropic transfer of energy 
to zonal flows observed in Chekhlov et al. [1] (See Section 8 for more discus- 
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0.04 
0.02 






Rh = ». 
Rh = 5 
Rh = 2 
Rh.1 
Rh = 0.5 



Fig. 8. Maximum growth rate max5R(cr(k; ap, i?e, i?/t)) vs ap at Re = 10 and 

|k|<l 

Rh = oo, 5, 2, 1, and 0.5. 



:a) Rh = oc 



(b) Rh 



(c) Rh 






Fig. 9. Contour plot of 3fi((7(k; ap, i?e, For each k, p is chosen to be 



aic - 7r/2. 



sion). The inflectional instability is no longer dominant for smaller Rhines 
number {Rh < 1/2 for example) when the resonant triad interaction becomes 
dominant. 

Figure 8 and 9 suggest that for an isotropic base flow (6) with large m, the 
strongest linear growth rate could be observed along = for large Rhines 
number {Rh > 2 for example). Thus, the modes with close to may play 
more important roles than the modes with ap. close to n/2. In Section 7, we 
consider the linear growth rate for more isotropic base flows (5) with m > 2 
by solving a severe truncation of (9) to verify this suggestion. 



5 Critical Reynolds number in Rh — > limit 



Here we consider the critical Reynolds number Re'^{ap, Rh) in the limit Rh 
0. The main result is 
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lim Re'^{ap, Rh) 

Rh — >Q 



V2, 



oo, 







< ftp < 7r/2 



(49) 



Q!r 



7r/2. 



For CKp = 7r/2, the critical Reynolds number is oo in the sense that for any 
Reynolds number, the base flow becomes stable if the Rhines number is small 
enough. Similarly, the critical Reynolds number is (or \/2) for < Op < 7r/2 
(or ftp = 0) in the sense that for any Reynolds number larger than (or 
\/2) the base flow is unstable if Rh is small enough. As discussed in Section 
4, resonant triad interactions are the dominant mechanisms for instability in 
the limit Rh 0. In fact, the result for < ctp < 7r/2 in (49) is obtained 
by considering wavevectors k resonant to p as shown later in this section. 
We will show that the expression for Re'^{ap, Rh) in the limit Rh — > is 
asymptotic to the critical Reynolds number obtained from a single resonant 
triad for < ctp < 7r/2. Thus, we flrst consider linear stability of a single triad 
interaction. 



5.1 A single resonant triad interaction and the critical Reynolds number 

Linear stability of a single triad k, p and — (k + p) with a forced wavevector 
p and a response wavevector k can be treated by a truncation of (12) with 
n — 0,1 (See [2]). This corresponds to 

ao + — = (50) 
obtained by a truncation of (15). Using (24), (50) can be written as 

2k{a + P/Re) k"^ sin a{2 cos a + k) 



sinQ;(l-A;2) 2q'^{a + q'^/Re) 



= 0, (51) 



where q = \/l + 2k cos a + k'^ = |k + p|. Taking the real part of (51), we flnd 
that the growth rate (7^ = 3?(f7) satisfles 



2/r>\/ 2/r>\ 2 (^r + /^^/-Re) (cT^ + g^/i?e) 



[ar + kyRe){ar + q /Re) + uf- 

z.Or -f- {K~ -f- q~) I ne 

(52) 



2(7^ + (A;2 + 52)/^e 
k sin^ q;(2 cos q; + /c) (1 — /c^) 



4g2 

The critical Reynolds number for instability is obtained by setting cr^ = in 
(52) and is denoted by Rej.\ Rel^ satisfles 
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_2 sin a(2 cosa + — ) uji 



f^""^]" " W ReUk' + q'Y ^^^^ 

If k is not resonant to p, nonzero Ui increases the critical Reynolds number. 
Thus, the smallest critical Reynolds number is observed at Rh = oo {ui = 0) 
for k not resonant to p. In the case that k is resonant to p, we have ui = 0, 
and (53) is reduced to 



rpc., p. ^2 4fc(l + 2fccosa + fc^)^ 

[i?e,(k;ap,m = oo)] = ,^,2 .(^ eosa + fc)(l - P) ^'^^ 
yielding the expression obtained in Smith and Waleffe [2]. 



5.2 < ttp < 7i/2 



For < Op < 7r/2, consider wavevectors k resonant to p and with ky > 0. We 
show (49) by deriving that, with a fixed Op, 



ro c/1 r,u ,sm2 Ak{l + 2k cos a + 

[Re'{k;ap,Rh = k')] = . / ^^77"^:^^ + ^ ^ 

sm Q;(2cosa + A;)(l — K^) 

4A;(l + 2fcsinap + A;2)2 
- ^ p ^ +o(/c), (56) 



cos^ a (2 sin ap + k){l — k"^ 



where s is a positive number depending on ctp to be determined later. The 
remainder terms o{k) depend on s and ttp. Relation (56) implies (49) since 
(56) approaches zero as A; ^ 0. The leading order term in (55) is the same 
expression for the critical Reynolds number (54) of a resonant triad, and (56) 
is immediate from (55) by (47). In Figure 10(a), the leading order term in (56) 
(solid lines) and numerically computed critical Reynolds numbers (symbols) 
are plotted for k = 0.1, 0.05, and 0.005 with s = 5 and < Op < tc/2. For 
each k, the lines and symbols are almost identical verifying (56). Figure 10(a) 
also verifies that the critical Reynolds numbers are decreased to for < Op < 
7r/2 as A; is decreased. In particular, the critical Reynolds numbers become less 
than a/2, the critical Reynolds number for the 2D Navier Stokes equations. 
An interesting observation is that the critical Reynolds number for Op = is 
2 in the limit /c — yielding a boundary layer at ap = 0. Also, for a fixed 
k, the critical Reynolds number in the limit ctp tt /2 is infinity resulting 
another boundary layer at Op = 7r/2. 



^ The Matlab function 'fzero' is applied to (12) to find critical Reynolds numbers 
numerically. 
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Figure 10(a) suggest that s = 5 is large enough to obtain (56) for < ttp < 
7r/2. In Figure 10(b), the critical Reynolds numbers are computed with s = 3, 
and agree with the leading order term in (56) except near ap = tt/Q. This 
observation can be explained as follows. From (23) and the fact that uji — 0, 
we have 



^-k^D^u{p)-k + 0{kyRh) 

= -(0, kyfD^u;{p) ■ (0, ky) + 0{k^/Rh). (57) 

The last equality is a consequence of k^ = 0{ky) along the resonant trace. 
One can show that for ctp = 7r/6, the first term in (57) is zero. Thus we have 



0(kyRh), ifap = 7r/6, 
w_i = <( (58) 

0{k^/Rh), otherwise. 

To have a;_i — > oo in the limit — > 0, we need s > 4 for ctp = tt/Q. For 
ttp 7^ 7r/6, s > 2 is enough to have ui^i oo. However, to derive (56), we 
need to assume that s > 4.5 for ctp — n/G and s > 2.5 otherwise, as we will see 
at the end of this section. For tt/G - 0.0001 < dip < n/6 + 0.004, Figure 10(c) 
shows the critical Reynolds numbers for k — 0.05 with s = 3, 4 and 5. The 
divergence near ap = 7r/6 is strong for s = 3, and almost disappears for s = 5. 
Relation (49) for < Op < it/2 is also verified in Figure 10(d) where the 
growth rates 3fJ((T(k; ctp, i?e, i?/i)) are obtained numerically from (12) along 
wavevectors k resonant to p with ap = 7t/9, tt/A and 7r/3. The horizontal line 
is the /Cy-coordinate of k. We choose a small Reynolds number Re — 0.1 and a 
small Rhines number Rh = 10~^, and a positive growth rate is observed along 
wavevectors k resonant to p for each ap. 

In the rest of this section, we sketch the derivation of (55). Since we are 
interested in /c <S 1, we use (18) and set 3?((T(k; ccp, i?e, Rh)) = to find the 
critical Reynolds number. Omitting details, we find that the critical Reynolds 
number Re'^{k; ap, Rh) satisfies 



iRe^'p + k 



2 [Re'']^ksm'^ a 



k + 2 cos a 



k — 2 cos a 



q'^{-l){ipRe^ + q'^{-l) + iu^iRe") 



g2(l)(ipi?e- + g2(i)) 

+ 0([i?ef A;^), 



(59) 



where p = '^{a{'k; ap, Re'^, Rh)) and the fact that k is resonant to p is used to 
set LOi — 0. Recall that 
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g^(±l) = |k±p|2 = 1± A;cosq; + A;^. 
By choosing s large enough in Rh = k^, suppose that 



a;_ii?e'^(k; ap, Rh = k^) 



^0{k'/^). (60) 



Then, the the second term in the bracket of the right-hand side of (59) becomes 
higher order than the first term. Thus, balancing the lowest order terms, we 
obtain 



2 [i^e'^j^/csin^ q;(/c + 2cosq;) , 2x 



iRe^p + k'^^ ir-—, ^ — +o(k^). (61) 

From the imaginary part of (61), one can see that Re'^p = o{k'^). Thus, the 
real part balance in (61) gives (55). Finally, we have Re%k; ctp, Rh) = 0{k^'^) 
from (55). Thus, by (58), (60) is true for s > 4.5 for ccp = 7r/6 and s > 2.5 
otherwise. 



5.3 dip = and — 7r/2 

For ctp = 0, in Figure 10(a) and (b), the critical Reynolds number along 
wavevectors resonant to p in the limit Rh — > is two, which larger than (49). 
In fact, we need to consider wavevectors of the form k = (0, ky). One can show 
that for sufficiently small \ky \ <^ 1, 

Re^iiO, ky); ap = 0, Rh) = V2. (62) 
for any fixed Rh. Relation (62) implies by definition that for any Rh, 

Re^ap = 0, Rh) < V2, 

which in turn imphes inequality "<" in (49). To show (62), note that (16) at 
k = (0, ky) is reduced to the growth rate for the 2D Navicr Stokes equations 
with k _L p since c^(k) = 0. Thus, the critical Reynolds number at k = (0, ky) 
is (see [9], [16] and [21]). It remains to show the other inequality ">" in 
(49) for ctp = 0. This requires rather complicated analysis, and is shown in 
Appendix A. 2. 

For ftp = tt/2, we know that the resonant triad interaction can not transfer 
energy from the wavevector p = (1,0) to other wavevectors ([7], [22], [23]). 
This fact can be also verified from (54). As shown in Figure 3, wavevectors k 
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(a) Re%k; ap, Rh = k') 




(c) Re%k; a^, Rh = k') 



?■ 5 
3 



O.S235 0.524 0.5245 0.525 0.5255 0.526 0.5265 0.527 0.5275 



Rh=k'' 
Rh=k'' 
Rh=k' 



(b) i?e<=(k; ap, Rh = k') 




(d) 3?((7(k; CKp, Re, Rh)) 





a .1/9 




— p 

_ _ . a -1/4 








- - p 








0.5 1 1 


5 i 


k 


x10^ 



Fig. 10. In all figures, k is resonant to p. In (a-c), the solid lines are the leading 
order critical Reynolds numbers in (56), and the symbols are numerically computed 
critical Reynolds numbers In (a) and (b) , the horizontal dotted lines are \/2. In (b) 
and (c), the vertical dotted lines are ap = 7r/6. (a) s = 5. (b) s = 3. (c) k = 0.05. 
In (d), the horizontal axis ky is the y-coordinate of k along resonant traces. 



resonant to p have ky 
obtain 



—1/2. Thus, we have cosa = —1/2. Prom (54), we 



[Re%k; Up = 7r/2, Rh 0)f 



Ak^ 



sm^a{l-k)^{k+l) 



< 



which is impossible. Thus, it is expected that the base flow becomes stable in 
the limit Rh when ctp = 7r/2 since the dominant instability mechanisms - 
resonant triad interactions - can not transfer energy to other wavevectors. We 
omit the derivation of (49) for ap = 7r/2, which is similar to that for ap = 
in Appendix A. 2. 
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6 Compcirison with Manfroi and Young's Results 

Manfroi and Young [20] studied the large-scale instability of a sinusoidal Kol- 
mogorov base flow ^'^s' given by 

^0^ = -cos(p-x). (63) 

With this base flow, the main result is that the critical Reynolds number in 
the limit Rh ^ oo is given by 

1^2, forap^O 

0, for < ftp < 7r/2 (64) 
"\/2, for Up — 77/2 

Here, i?e^y (k; ap, Rh) is the Reynolds number which makes the leading order 
term for the growth rate equal to zero. (The formula for the growth rate by 
Manfroi and Young [20] is given by (71), derived at the end of this section.) 
i?e^y (k; ap, Rh) is analogous to i?e°(k; ap, Rh) in Section 3 and is understood 
as the critical Reynolds number in the hmit A; — > 0. Note that Manfroi and 
Young use the nondimensional number which is same as Re/Rh in this 
paper. 

For < CKp < 7r/2 in (64), an interesting observation is that the critical 
Reynolds in the limit Rh — cxd is less than y/2, the critical Reynolds number 
for the case Rh = oo. In other words, there is a discontinuity in the critical 
Reynolds number at Rh = oo. The result (64) is different from our result for 
a Rossby wave base flow ^/t(x, t; p) given by 

^'?j(x,t;p) = -cos(p-x-a;(p)t). (65) 

For (65), it is observed that the marginal Reynolds number Re^{]s.; ap, Rh) 
satisfies (33) for any ap. Thus, there is no discontinuity in the critical Reynolds 
number in the limit Rh — > oo for (65). 

To see the difference between two base flows ^/j) consider the forces 

F/^ and Fr that respectively maintain the base flows. For the Kolmogorov 
base flow ^^C) the force F^ satisfies 

l-(V.F,)..^^{^)\{l)\4p-^-x) 



lim mini?e^y(k; CKp, i?/i) = < 

Rh—* oo k 



34 



and X = tan ^(pxRe/Rh) (66) 
For the Rossby wave base flow the force Fr satisfies 

i-(V X F«) • z = i- cos(p • X - u;{p)t). (67) 

From (66) the phases between F^: and different by x if p-^ 7^ {ap 7^ 

7r/2) and < 00. In particular, since (66) does not hold for i?e = 00 ( x 
becomes infinite at Re — 00), ^'^c ^^t a solution of the unforced, inviscid 
/3-plane equation. In contrast, the force Fr maintaining the Rossby wave base 
flow \E'5j has the same phase as ^jj, and is a solution of the unforced, inviscid 
/3-plane equation. In the case ctp — 7r/2, two base flows (63) and (65) are 
identical and the critical Reynolds number is y/2 for both cases in the limit 
Rh — > 00. 

The discontinuity in the critical Reynolds number at Rh = 00 for base flow \l'5f 
is probably linked to the fact that "if^ is not a solution to the inviscid, unforced 
/9-plane equation. Manfroi and Young pointed out that a discontinuity in the 
critical Reynolds number is observed in the 2D Navier-Stokes equations with 
base flow 



^%{x.;6) — sin X + cos{Sy) . 

In [26], Gotoh and Yamada showed that the critical Reynolds number for 6=1 
is different from the critical Reynolds number in the limit 5 — > 1. Similar to 
^5c) for (5 7^ 1 the base flow is not a solution of the 2D unforced Euler 
equation. This is immediate since J{V'^^%, 7^ for (5 7^ 1. 

In the rest of this section, we sketch the derivation of large-scale growth rates 
for both ^^^f ^iid by Multiple Scales analysis following Manfroi and Young 
[20]. Further details are given in Appendix A. 3. To this end, we consider a 
base flow 



*m(x, t; p, c) = - cos(p • X - cuj{p)t). (68) 

Note that t; p, 0) = *?f(x; p) and t; p, 1) = *?j(x, t; p). Follow- 

ing Manfroi and Young, we consider a rotation of the coordinates by ap. As a 
result, the base flow becomes periodic in x, and the Multiple Scales analysis 
becomes simpler. 

With rotation by ctp, the /3-plane equation is given by 



(V2^)i + J(VV,^) + il/Rh,)^l^, - {l/Rhy)i^, 
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(l/i?e)VV + (l/^e)(V X F) • z 



(69) 



and the base flow becomes 



K 



(x, t] p, c) 



cos(x — cc<;((l, 0); a-p)t) 



(70) 



where {1/Rhx,l/Rhy) — (l/i?/i)(cos ctp, sin ctp). The rotated dispersion rela- 
tion a;(k; ctp) is 



We consider scahngs similar to [20]: 

dt + edr + e'^dr, d^^d-^ + edx, dy ^ edy 
Rh = e-^Rh\ Rh^ = e'^Rhl, Rhy = e'^Rhl 

-0 = ^0 + + £2^2 _^ _ 

Our interest is in large Rhines number Rh — e~^Rh^, and large scales X = 
(X, Y). At such spatial scales, the growth rate is expected to be 0{e^), thus the 
O(e^) time scale r is introduced (see [6], [16], [14], [20], and [15]). Large scales 
correspond to wavevectors k with k = 0{e). Thus, at large scales the disper- 
sion relation, which is an inverse time scale, satisfies a;(k; ap) = {—Rhxkx + 
Rhyky)/k'^ — 0{1) since Rh — e'^Rh^. This motivates the 0(1) time scale t. 
Finally, since the base flow depends on the time scale Rh~H = 0{s), the 
0{e) time scale T is introduced. Since the Kolmogorov flow (63) does not 
depend on T, it is sufficient to consider dt —>■ dt + e^d-j- in Manfroi and Young 



a;(k; ccp) 



[20] 



To find the linear growth rate, we set the leading order term ijj^ to 




where the large-scale dispersion relation is 



Kx/Rhj + Ky/Rhl 
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Applying Multiple Scales analysis, we find that the coefficient A{t) satisfies 



At^O and = aA, 
for a complex number a. The real part of a is the growth rate and is given by 

3f?((7(k; ttp, Re, Rh} , c)) = -— (1 

Re 

Kl - TKl + {KVlRef - 4K^nRe\-c + l)/Rhl 

where fl = fl(K;ap). The derivation of this growth rate is very similar to 
the one found in Manfroi and Young [20], and the main steps are given in 
Appendix A. 3. 

The expression (71) with c = is Manfroi and Young's result for the large-scale 
growth rate. Then, i?e^jy (k; ap, Rh) in (64) is found by setting the growth rate 
(71) with c = equal to zero. By taking a special path for k — s> and Rh — > 
for each ap, Manfroi and Young [20] obtained the critical Reynolds number 
(64). For c — 1, after adjusting the scales and rotating the expression by ap, 
(71) is identical with the leading order term in the growth rate (16) obtained 
from the continued fraction. Mathematically, it is the term 4Kx^Re'^{—c + 
1)1 Rh\ that leads to the difference between growth rates for base flows ^'jc 
and 

In earher work of Manfroi and Young [15], they considered a base flow of the 
form 




= COSX+ C/y. (72) 

This base flow is obtained by adding a term Uy to the Kolmogorov base flow 
(63) with ftp = 0. As Manfroi and Young [15] pointed out, this term accounts 
for a Galilean translation in the x-dircction with the speed U. Through this 
Galilean translation the Rossby wave base flow (65) becomes (72) if 1/ Rh = U . 
Under this translation, one can show that the /9-plane equation remains same. 
Thus, the base flows (72) and (65) are equivalent under the translation if 
1/Rh — U. With base flow (72), Manfroi and Young obtained the critical 
Reynolds number a/2 if 1/Rh = U. This is consistent with our result that the 
critical Reynolds number is y/2 with the Rossby wave base flow (65). 



7 Bcise flows with more than one Rossby wave 



37 



7.1 Growth rates and critical Reynolds numbers for the special flows with 
m = 2 or m = 3 

In this section, we consider large-scale (fc -C 1) growth rates for two special 
base flows whose stream functions are given by 

2 

^'2(x,t;pi,p2) = -^cos6'(x,t;pj) where Op^ - Op^ = 7r/2, (73) 
i=i 

^3(x,t;Pi,P2,P3) = -Y.cose{x,t;pj) where ap^^, - = — . (74) 

Note that the base flow (73) has two Fourier wavevectors pi and p2 (and their 
conjugates) and the base flow (74) has three wavevectors pi, p2 and ps (and 
their conjugates). The choice of these base flows is motivated by the work by 
Sivashinsky and Yakhot [16] of the large-scale instability of the following base 
flows in the 2D Navier Stokes equations 



$jL = cosx 

$2 = cosxcosy = -(cos(x + y) + cos{x — y)) 

$3 = sin 2x — 2 sin X sin ^y = sin(2x) + cos(x + VSy) — cos(x — \/3y). 

Similar to (73) and (74), $2 ^^^1 $3 have two and three Fourier wavevectors 
(and their conjugates), respectively. Using Multiple Scales analysis, Sivashin- 
sky and Yakhot showed that base flows ^2 have the critical Reynolds 
number v^-[f] An interesting result is that for the base flow $3, the linear 
growth rate at k with <C 1 is independent of the direction of k (i.e., is 
isotropic) and is negative for all Reynolds numbers. Thus, there is no direct 
instability to large-scales for the flow <I>3. In this section, the same results is 
obtained for the flow (74) with Rh = 00. 

The main result of this section is a generalization of [16] to the /5-plane equa- 
tion for the base flows (73) and (74). Sivashinsky and Yakhot obtain the 
large-scale linear growth rate using Multiple Scales analysis with the follow- 
ing scaling: 



The definition of the Reynolds number for ^2 Sivashinsky and Yakhot is dif- 
ferent from our definition. In their definition, the critical Reynolds number for <I>2 
is 2V2. 
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9x = + edx, dy^dy + edy, dt = e'^dr- 



Since the base flows $2 ^^^1 $3 depend on both x and y, the analysis becomes 
rather complicated. Later, Waleffe [21] obtained the same result easily by 
direct apphcation of formula (16) for m = 1 with Rh — 00. We follow Waleffe's 
idea to find the large-scale growth rate for the base flows (73) and (74), and 
verify the results by Multiple Scales analysis in Appendix A. 3. 

Let (73) be the base flow which contains two Fourier wavevectors (and their 
conjugates). Without loss of generality, we can assume that < ctp^ < 7r/4. 
For each wavevector k, let aj = — ftp^ for j = 1, 2. Since apj— api = 7r/2. wc 
have 0:2 = ai — 7r/2. By adding terms directly from (16) with p = Pj, j = 1, 2, 
the growth rate is given by 



3?((T(k; ttpj, ap2, Re, Rh)) 
Re^ 



Re 
Re^ 



1 



sm «! 



2 {l + Re^Lo^y 



r{l-8cos^ai + i?eV} 



sin^ a2 



2 (l + i?e2a;2)2 



{1 -8cos^ a2 + ReW} 



+ 0{k^) 



Re 



1 „ 2I + ^e2a;2 -4sin^(2ai) 
1 — Re - 



2{1 + Re^uj^y 



+ 0{k^). 



(75) 



Note that we only consider the growth rate up to the order of 0{k'^). This is 
Waleffe's method [21]. In particular, for Rh = 00 (a; = 0), (75) becomes 



3ft(cr(k; CKpj, ctpj. Re, Rh — 00)) 



Re 



I -Re' 



4sin2(2ai) 



+ 0(A;^). 



(76) 



The higher order term is 0{k^) since the third order term in (16) is zero for 
Rh = 00. This expression (76) for large-scale growth rate was obtained by 
Sivashinsky and Yakhot [16] and by Waleffe [21]. 

Similar to Section 3.2, for base flow (73) we define the marginal Reynolds 
number Re^ (k; Op-^^ , , Rh) by the smallest nonnegativc Reynolds number 
which makes the leading order term of the growth rate (75) equal to zero by 
solving 



2{1 + Re^uj^)^ 



Re' 



l + i?e2a;2-4sin^2Q;i 



and 1 + ReW - 4 sin^ 2ai > 0. (77) 
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If there exists no Reynolds number satisfying (77), the marginal Reynolds 
number is defined to be infinity. i?e°(k; ctp^, ctp^, -R/i) is understood as the 
critical Reynolds number in the limit k ^ 0. For Rh = 00 {u = 0), from (77), 
the marginal Reynolds number is given by 



i?e°(k; ttpj, Opj, = 00) = w 5 , if 1 — 4 sin^ 2ai > 0, 

V 1 — 4 sin 2ai 

and its minimum is a/2 when ai = or 7r/2. Thus, if Rh = 00, the critical 
Reynolds number is \/2 for the 2D Navier Stokes equations with base flow 
(73). This is the same result of Sivashinsky and Yakhot [16] and Waleffe [21] 
As in Section 3.2, we see that i?e°(k; Op^, ap^, Rh) is always larger than or 
equal to \/2, since we have 



Re {k;ap,,ap^,Rh} -2 = .,-022 — a ■ 20 0' 

if (1 + Re^uj"^ — 4sin^ 2ai) > 0. For ai = or 7r/2, we also have 



lim i?e°(k;ap,,ap2,m) = 72, (78) 

recovering the critical Reynolds number a/2 for the 2D Navier Stokes equations 
in the limit Rh —>■ 00. 

Similar to Section 3.3, for < ap^ < 7r/4 we find that the critical Reynolds 
number in the limit A; — is given by 



lim min Re (k; Up-^, ap^, Rh) 

ko-*0 |k|<fco 

\ t7 T < an, < ■^sin~^(— 7=) 

V 1 -4sin2(2apJ - pi - 2 

(79) 

4v^|sin(2apJ| I sm~\^^) < ap, < ^ 

We omit the derivation since it is much like that of (34) for m = 1. Similar to 
(34), (79) does not depend on Rh. 

Now, consider the base flow (74) which contains three Fourier wavevectors 
(and their conjugates). For each wavevector, let aj = — ap. for j = 1, 2, 3. 

^ See the footnote of 4. 
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Since a2 = cci — 27r/3 and 0:3 = ai — 47r/3, we have two simple trigonometric 
identities: 

3 

^sin^ aj = 3/2 

3 

sin^ aj cos^ aj — 3/8. 

Using these identities, by adding terms from (16) for ccp = ccp^., j — 1,2,3, 
the growth rate is given by 



K(cr(k; ctp^, ctpj, ctpg, i?e, Rh)) 



Re 



Re^ 



2(1 + Re^u'^y 



;i + i?eV)^sin2 



-8 ^ sin^ Qij cos^ aj 



+ 0(A;=^) 



Re 



3Re' Re^u^ 



4 (1 + Re^uj^)^ 
In particular, for = 00 it becomes 



+ 0{k^). 



(80) 



3?((7(k; CKp^, CKpj, Q;p3, i?e, = 00)) = — — ( 1 + -Re^ ] + 0{k^) 

Re V 4 / 

and the leading order expression is isotropic in k and always negative. Thus, 
there is no direct instability to large scales for all Reynolds numbers. This is 
the result by Sivashinsky and Yakhot for the "isotropic" base flow $3 and also 
obtained by Waleffe [21]. 

An interesting observation is that for Rh < 00, the large-scale growth rate 
(80) does not depend on the particular orientation of pj. However, (80) is not 
isotropic in k because the /3-plane equation itself is not isotropic. The growth 
rate (80) depends on k only through a;(k), whereas growth rates (16) and (75) 
depend on k not only through a;(k) but also through a or cci, respectively. 

As in Section 3.2 and above for (73), we define the marginal Reynolds num- 
ber i?e° (k; ttp^ , ttpj , Q;p3 , Rh) by the smallest nonnegative Reynolds number 
satisfying 
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R^" = l ^^ + ^f"^7 ^ and Re'u' - 1 > 0. (81) 
3 Re^cu^ — 1 

When Rh = oo (cu = 0), there is no nonnegative Reynolds number satisfying 
(81). Thus, the marginal Reynolds number is infinite and there is no direct 
large-scale instability [16], [21]. Analogous to (34) and (79), one can show that 



niini?e°(k; Up^, apr^, ccpg, Rh) — ■\/32/3 (82) 

for Rh < oo. We omit the derivation of (82) since it is similar to that of (34) 
for m = 1. The derivation uses the fact that 



min i?e°(k; ttp^, Q;p2, Q;p3, i?/i) = min Re^ (k; ap^^, ap^, ap^, Rh) 

k^O |k|<A;o 

for any positive number ko, since (81) does not depend on cti. Another ingre- 
dient in the derivation of (82) is that the function 



m = lj^^ for x>0 (83) 
has maximum value 3/32 at x = 3. We conclude that there is a direct instabil- 



ity to large scales (A; <C 1) with the "isotropic" base flow (74) for Re > ^32/3 
and Rh < oo. This is different from the 2D Navier Stokes equations, where no 
direct large-scale instability is observed for any Reynolds number with base 
flow (74) [16]. 



7.2 Numerical study of the growth rate for various Rh 



In this section, we consider the linear growth rate at wavevector k with k < 1 
for large-scale instability with m > 2 in (6). As discussed in Section 2, it is not 
straightforward to solve (9) numerically for m > 2 since the technique of the 
continued fractions used for m = 1 can not be generalized. However, in Section 
4.2 we also observe that for m = 1 the severe truncation (48) with n — —1, 0, 1 
is a qualitatively good approximation for small k. With this background, we 
truncate (9) to a finite dimensional eigenvalue problem with 

n = 0, ±ej, for J = 1, • • • , m. 
In terms of triads, this truncation involves only 2m-triads: 

(k,Pj,-(k + Pj)), (k, -pj,-(k-pj)), j = l,---,m. 
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For example, for m = 2 the corresponding eigenvalue problem with the trun- 
cation (9) is given by 



-c(0,-l) 



62(0,0) 



-ci-l,0) ^ 



fa2(0,-l) fel(-l,0) ((. bl(l,0) 62(0,1) 

■ g2(0,0) g2(o,o) ^-V^' g2(o,o) (/2(0,0)2 



fei(0,0) 
g2(l,0) 



b2(0,0) 
■92(0,1) 



-c(l,0) 



-c(0,l) 



where 



c(n) = ;^g2^n) +iu;„, 



&.(n) = -2(P. xq(n))(g2(n)-l), 



, * = (V'(-1,0), -0(0,-1), -0(0,0), -0(1,0), '0(0,1))^- 

Moreover, we consider for each m, the "isotropic" base flow with 



«P,=J — , J = !,•••, ^■ 
m 

For large m, this base flow mimics deterministic forcing over a shell of radius 
one in Fourier space. For each wavevector k with Re and Rh fixed, we find the 
eigenvalue of the truncated eigenvalue problem with the largest real part, de- 
noted by a(k; m, Re, Rh). By the "isotropic" choice of cip^. , 3?(a(k; m. Re, Rh)) 
is symmetric about both the x-axis and the y-axis if m is even. 

Each contour plot is for wavevectors k with 0</Ca;<l,0<A;j/<l and |k| < 1 
to study the large-scale instability. We use a grid of size 100 x 100 over < 
kx < 1 and < ky < 1. Figure 11 shows the growth rate ^{a(k;m. Re, Rh)) 
for m = 30 and Re = 10. The growth rate for Rh = cxo is invariant under 
the rotation of 7r/m since the Navier Stokes equations are isotropic. Also, note 
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{a) m = 30, Rh = oo {h) m = 30, Rh = 2 (c) m = 30, Rh = 1/2 




Fig. 11. 3f?((7(k; m. Re, Rh)) at m = 30 and Re = 10. 3f?((7(k; m. Re, Rh)) is computed 
from (9) with the severe truncation. 



(a) m = 30, Re = 1000, Rh = 2 (b) m = 30, Re = 1000, Rh = 1/2 




k 



Fig. 12. 3fi(a(k; m. Re, Rh)) at m = 30 with Re = 1000 and m = 40 with i?e = 10. 

that there is no direct instabihty into large scales, as consistent with the result 
in Section 7.1 for the "isotropic" flow (74) with m — 3. For Rhines number 

2 < Rh < oo, the largest growth rate is observed to be exactly along the line 
kx = 0. However, for Rh = 1/2, the largest growth rate shifts slightly away 
from kx = 0. Again, this transition can be understood in terms of resonant 
triad interactions, which cannot transfer energy into wavevectors with kx — 0. 
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A similar result transition is observed in Section 4.3. The 0(1) Rhines number 
for the transition seems to be independent of Re and m, as demonstrated in 
Figure 12 for Re = 1000 with m = 30, and Re = 10 with m = 40. 



8 Summciry 

The growth rate of linear instability at large scales is obtained for Rossby 
wave base flows maintained by deterministic forcing of one or several Fourier 
modes. For forcing of a single mode, we use the technique of continued fractions 
following [9], [19], [10], [11] and [12], and asymptotic analysis. We show that the 
mechanism for linear instability changes from inflectional to triad resonance 
at an transition Rhines number, independent of the Reynolds number. 
The critical Reynolds number for instability has also been obtained in various 
limits. In the limits k ^ and Rh oo, we recover the classical result 
Re'^ = \/2 found by Meshalkin and Sinai [9] for the 2D Navier Stokes equations 
(see also [16], [14], [15] and [21]). Furthermore, we generalize [9] to find Re'^ 
given by equation (34) for A; —> and finite Rhines number Rh. For Rh — > 0, 
we find Re"^ = for Rossby wave base fiows of all orientations except zonal 
and meridional {i.e., for < < 7r/2, where ap measures the angle between 
the forcing wavevector p and the zonal direction x, and |p| = 1). For Rh 0, 
a zonal base fiow is stable and a meridional base fiow has Re'^ — V^. It is clear 
that, for small Rhines numbers, resonant triad interactions reduce the critical 
Reynolds number below -\/2 for < ctp < 7r/2. 

In order to study more isotropic forcings, we consider base flows consisting of 
m > 1 Rossby waves, thereby generalizing [16] to the /3-plane with Rh < oo. 
For m = 3, we flnd instabihty for Re > ^32/3, whereas the analogous base 
flow was found to be stable for Rh = oo in [16]. For large m, such a base 
flow mimics deterministic forcing over a shell of radius one in Fourier space. 
Numerical computations of growth rates for m = 30 and m = 40 show that the 
most unstable mode is purely zonal with /c^, = for 2 < Rh -C oo (Figure 11 
and 12). This is because one of the forced Rossby waves is meridional with 
CKp = 0, and the inflectional instability is strongest for meridional base flows, 
as shown in Figure 8 and 9, and discussed in Section 3.3. As the Rhines 
number is decreased from Rh = 2 to Rh = 1/2, the most unstable mode k 
shifts to a nearly zonal flow with not exactly zero (Figure 11 and 12). This 
transition in the most unstable wavevector can be understood as the change 
from inflectional instability to resonant triad instability at Rh = 0{1), and 
reflects the result that resonant triads cannot transfer energy directly to zonal 
flows with k^ = [1] (see also [7], [8], [6], [2], [5] and Section 3). Waleffe 
[22] showed analytically that resonant triad interactions can transfer energy 
to nearly zonal flows with near zero. In Figure 11 and 12, however, one 
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also sees that wavevectors with kx exactly equal to zero remain unstable even 
for Rh = 1/2, though they are not the most unstable wavevectors. This is 
because for Rh = 0(1), the mechanism for instability is a combination of 
both the inflectional and triad resonance instabilities. The latter observation 
may be relevant to recent simulations of /3-plane flow driven by isotropic, 
stochastic forcing in a wavenumber shell [1], [2], [3]. With energy input rate 
€f and peak wavenumber A:/ of the stochastic force, the Rhines number for 
these simulations is appropriately defined as Rh = {efkjy^^ / P, measuring the 
relative strength of the forcing and P terms. In [1] , this definition of Rh leads 
to three values Rh — oo, Rh — 5 and Rh = 1. In the latter two cases, energy 
is transferred anisotropically to large-scale zonal flows with kx = 0. Even 
though our Rhines number is defined differently, it is nevertheless intriguing 
that linear instability mechanisms transfer energy to zonal flows with k^ = 
for our Rh — 2. Indeed, for Rh — 2, the most unstable modes are purely zonal 
flows. Thus it is possible that linear instability plays a role in the population 
of zonal flows observed in numerical simulations of the (nonlinear) /5-plane 
system at moderate Rhines numbers Rh > 1. Further study of stability of the 
stochastically forced system and of nonlinear transfer mechanisms is necessary 
to fully understand numerical simulations at both large and small values of 
the Rhines number. Finally we note that, for the atmosphere and the oceans, 
Rhines number 0(1) is associated with the most energetic eddies [19]. From 
the peak wavenumber of the atmospheric spectrum. Gill [19] estimates that 
the Rhines number of the atmosphere at 40''N and 500 mb is approximately 
0.9. 

Acknowledgements The authors thank Fabian Waleffe for the suggestion to 
use the method of continued fractions and for many thoughtful discussions. 
The support of NSF is gratefully acknowledged, under grant DMS-0071937. 



A Appendix 



A.l On the continued fraction 



We briefly discuss the convergence of the continued fractions in equation (15). 
To this end, for complex numbers ^„ and rjn, consider the continued fraction 



K=—^^. (A.l) 



mA ^ 



We have the following classical theorem on the continued fraction by Sleszyhski 
and Pringsheim [27]. 
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Theorem The continued fraction (A.l) converges if for all n 



\Vn\ > l^nl + 1- 

Under the same condition, we have 
\K\ < 1. 



In equation (15), we have two continued fractions with = 1 and rjn = a±n. 
Since we consider the case Re < oo, for fixed a, k, ctp, Re and Rh we have 

hm |r7n| = hm \a±n\ = oo. 

n— >oo n— >oo 

Thus, by the theorem above, both continued fractions converge. Next, we show 
(21) for A; <C 1. Since we have 



111 

< CRe k<- for A; < — — and n = ±2, ±3, . . . , 



|a„| - - 2 - 2CRe 

from (20), by the theorem, we have 



Q±3 + a^,+^ 



< 1. (A.2) 



Prom (A.2), one can derive (21) using an asymptotic expansion in k and (22) 
is a consequence of (21). 



A.2 Critical Reynolds numbers for ap — and a-p — t: /2 in the limit Rh — > 



In this appendix, we consider (49) for ctp = and = 7i/2. Since (49) can be 
shown in a similar way for both cases, we show (49) only in the case ap — 
and omit the other case. 

Consider (49) for the case ctp = 0. Since we showed that 



lim Re^iap = 0, Rh) < ^2. 

Rh — >0 

in Section 5.3, it is enough to show that 
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lim Re'^iap = 0, Rh) > V2. (A.3) 

Rh^O 

For this end, we assume that for some k° we have 



hminf Re%k; = 0, Rh) = lim Re^a^ = 0, Rh) (A.4) 

By the symmetries of 9ft(cr(k; ctp = 0, Re, Rh)) discussed in Section 2, we can 
assume that < < 1/2 and ky > 0. We show (A.4) separately in each of 
two cases: k° = and k'^ ^ 0. 

First, consider the case k° ^ 0. If k*^ is not resonant to ±p, then we have 
\uj±i\ — >• cx) in the k — > k° and Rh hmits. Thus, (15) is reduced to 



ao = 
which gives that 



hminf a(k; = 0,Re,Rh) = -ik'y/Re < 0. 
This in turn imphes that 



hminf i?e^(k: = 0, Rh) = oo. 

If k° is resonant to — p = (—1, 0), then k° is not resonant to p = (1, 0). Thus, 
in the Rh — ^ and k ^ k° limits, we have uj-i —>■ uj^i for some uj^i, and 
|a;i| ^ cxD. Then, (15) is reduced to 



1 

ao H = 0. 

From this we can find the critical Reynolds number 



liminf \Re^(k] a„ = 0, Rh)f 



4A;0(1 -2A;0 costtk" + 
sin^ ako (-2 cos + A;0) (1 - ( fco) 



by (53). We claim that the right-hand side of (A. 5) is larger than or equal to 
4. This can be shown by closer analysis of (45) for Op = and we omit the 
proof. This completes the demonstration of (A.3) when k° ^ 0. 

^ Since < < 1/2, k'^ can not be resonant to p = (1,0). 
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Now, we consider the case k° — 0. Consider the following four cases of the 
limits k ^ and Rh 0. Here, e and 5 are real constants with for < £ < 
d<l. 

(1) < 2k^+' , 

(2) k^ > A;2+72 and k/Rh < 2k^ , 

(3) k^ > A;2+72, k/Rh > k^ and k^ < ky/V^ , 

(4) k^>ky. 

Note that these four cases are sufficient to cover all limits k — > and Rh — > 
in (A. 4). Of course, this is a formal argument to deal with the limits k — and 
Rh — > 0. To be mathematically rigorous, these limits have to be understood 
in terms of sequences. However, using sequences complicates the notation. 

Since we are interested in /c -C 1, we can use (18) and set 3fJ(cr(k; ctp. Re, Rh)) — 
to get 



,9,„^ A;sin^a(A; + 2cosa 

ip + k^/Re" - ^ 



4g2(l)(ip + g2(l)/i?e^ + io;i) 

k sin^ a{k — 2 cos a) 
Aq^{-l){ip+q^{-l)/Re'' + iuj_i) 



Oik"") (A.6) 



where p = Q{o-(k; ap, Re, Rh)) is the imaginary part of a{k; ap. Re, Rh), and 
g(±l) = |q(±l)| = |k±p|. Here, Re^ = i?e^(k; Op = 0,Rh). Prom (A.6), one 
can see that p(k) = 0{k) for A; 1 and all Rh. 

Now, we consider each case of the k — > and Rh — > hmits. 

(1) k, < 2k''+\ 

k 2k^'^^ 

Since cos a = < — ; — = 2k^~^^, the real part of (A.6) becomes 
k k 



e/Re^^k^^H 



sin^ a 



Aq^{l){ip + q'^{l)/Re'' + iuJi) 



^iq^{-l){tp + q^{-l)/Re'^ + zuj.i)J^ ^ ^ 

- " \Aq{l)^/Re- + Aq{-1)^ / Re'^ ) ^ 
Thus, in the k — > and Rh — > hmits, we have 
liminf Re%k; ap = 0, Rh) > V2. 

(2) k^ > A;2+72 and k/Rh < 2k^. 
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Using an asymptotic expansion of uj±i in and ky ior k <ti 1, we can 
show that 



^±^ = -^[Y-+k^±{k^-2kl)+0{k') 



(A.7) 
(A.8) 



We set 7(k) = k^/k"^ + k^. Then, similar to the derivation of (16), we can 
show that 



3fJ(cT(k; ctp. Re, Rh)) = — ^A;^ 

Re 



1 



2 sin^ a 



I Up- 

2 (1 + 72i?e2)2 



X 



{l-8cos'a + 7W}] +0(A;'+'^). (A.9) 

This is the same as (16) at leading order. In the k — > and Rh 
limits, the critical Reynolds number is obtained by setting the leading 
order term equal to zero. Similar to (31), we have 

liminf ReH\i; = 0, Rh) > V2. 

(3) k^ > A;2+72, k/Rh > k^/2 and k^ < ky/V^ . 
bmce — > -—TTT = — , we have 



A;2 - 2 A:2 
k±i(k)| 



- Rh" 



where Ci is a constant independent of Rh. Combining this with k/Rh > 
{l/2)k\ we have 

k±i(k)|>^F>(Ci/2)A;-^+^+^ 

Note that since < ky/-\/2>, we have k"^ — 2kl > k'^/2. Using this and 
(A.7), one can show that 

|<^i| > for k <^1. 

Since we have |p(k)| = 0{k) and |a;±i| > {Ci/2)k~^~^^~^'^ , (A. 6) becomes 

k sin^ a{k + 2 cos a) 



ip + k'/Re' 



Aq^{l){q^{l)/Re- + iuJi)[l + 



ip 



q^{l)/Re^ + iu;i^ 
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+ 



k sin^ a{k — 2 cos a) 



V(-l)(g2(-l)/i?e- + ia;_,) 1 + 



'2(-l)/i?e^ + iu;_i 



) 



k sin^ a cos a 



2q''{l){q''{l)/Re- + iu;i) 





(A.IO) 



2?2(-l)(g2(-l)/i?e- + a;-i) 



We note > q{—l) and |cc;i| > |ti;-i|. Taking the real part from each 
side of (A.IO), in the /c — > hmit we have 

hmsup 1/Re''<0. 

This in turn imphes 

hminf Re^ik; = 0, Rh) > V2. 



One can see that 

C 

\u)±i\ > — rk~^ for some constant Co 
' ' - Rh 

Then, from (A. 6), taking the real parts, 

k'^/Re" < CsRh k'^ + 0{k^) for some constant Cg 
In the k — > and Rh — > hmits, 

liminf i?e^(k; = 0, Rh) > v^. 

This completes (A. 3) when k° = and shows (49) for — 0. 
A. 3 Multiple scale analysis 

In this section, we consider the base flows (68), (73) and (74) and find the 
large scale growth rate of perturbations using Multiple Scales analysis. These 
base flows are spatially periodic. For Multiple Scales analysis, we rotate these 
base flows to be periodic in both x and y. Being periodic in the coordinate 
directions simplifies Multiple Scales analysis. For (73) and (74), let p = Pi- 
With the rotation of ttp, the /5-plane equation becomes (69) and the base flows 
(68), (73) and (74) becomes 



(4) k^ > ky/2 
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= -cos(p-x-ca;(p;Q;p)i) p = (1, 0) (A.ll) 

2 

*° = -^cos(p,-x-a;(p,;ap)t) pi = (l,0), P2 = (0, 1) (A.12) 

*3 = - cos(pj • X - uj{pf, ap)t) pj = (cos((j - 1) — ), 

sm((j - l)f )). (A.13) 

These are the base flows studied by Sivashinsky and Yakhot [16] for the large 

scale instability in the two dimensional Navicr-Stokcs equation {Rh = oo). 
Also, Manfroi and Young [20] considered the base flow (A.ll) with c — 1. 
Note that (A.ll), (A.12) and (A.13) satisfy 

V^^o + *o = 0, 
and with c—1 they satisfy 

(V2^°)t + (l/i?/i.)*° - {l/Rhy)^l = 0. (A.14) 

Recall that {1/Rhx,l/Rhy) — {1/ Rh) {cos Up, sin ap). A small perturbation 
^ = ^ - ^0 satisfles 

( VV)t + *;(V^V' + - *°(V^V + 

+{l/Rhx)iPx - il/Rhy)iPy = (l/i?e)VV (A.15) 

Following Sivashinsky and Yakhot [16] and Manfroi and Young [20], we intro- 
duce the scales 

dt^ dt + edr + e^dr, 8x^8^ + edx, dy^dy + edy, 
Rhx — s ^Rh^ Rhy — £ ^Rhy. 

The choice of these scales is discussed in Section 6. With these scales, (A.15) 
becomes 

{dt + edr + e'a,) + ^^[(V'V' + + ^(VV + ^)x] 

-^°[(V^^ + ^l^)x + eiV'^l^ + ^)y] + £/5^(^, + e^x) 

-ef^litl^y + eiI^y) = (l/i?e)VV (A.16) 
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xxxa; ~l~ "^^xxyy ~l~ ^yyyy) ~l~ 4£('0j;j;j;X + '4'xxyY ~l~ ^xyyX ~l~ ^yyyY^ 



+2e'^{3lp^:,XX + SV'j^yyy + V'a:a;yy + V^j/yXX + "^IpxyXv) 



+4£^(^a;xxx + V'xxyy + V'yxxy + V'^yyy) 



+ 2'^xxYY + V'yyyy)- 



Integrating (A. 16) in x and |/ over the periodic domain of we obtain 



Note that from (A. 16) and (A.17), *° depends on X, y, Rh^ and T. For the 
details of derivations of (A. 16) and (A.17), we refer Sivashinsky and Yakhot 
[16] and Manfroi and Young [20] (see also [6], [14] and [15]). 

We seek the solution of (A. 16) of the form 

-0 = ^0 + + ^2^2 _^ _ 

and for simplicity we assume that the base flow satisfies (A. 14). This assump- 
tion is not true for the base flow (A. 11) if c 7^ 1. 

Prom the zeroth (0(1)) approximations of (A. 16) and (A.17), we assume that 

V'° = ^(T,T)e^(^-^-^(^'"-)*) 
where the dispersion relation is given by 








(A.17) 



-RhjK^ + RhlKy 
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In the first {0{e)) approximation of (A. 16), we find 



Prom the first {0(e)) approximation of (A. 17), we have At — and thus 
A — A{t) in ■0°. Prom the second (0(£^)) approximation, neglecting small 
terms of the order 0{Re'^), we find 



Neglecting terms of the order 0{Re^) in the second {0{e'^)) approximation is 
discussed in Sivashinsky and Yakhot[16]. Then, the second approximation in 
the integral relation (A. 17) gives 



+dy J J ^S[2(€x + ^^y ) + i^xx + rYY)]dxdy 



(A.18) 



+Re ^ J J il^xxxx + SV'xxyy + ^I^yyyy 
Inserting ipQ, ipi and ip2, it reduces to 

At- = a A for some a. 

and the linear growth rate is obtained by taking the real part of a. We refer 
to Sivashinsky and Yakhot [16] for the details. 

In the case of the base flow (A. 12), (A.18) gives the growth rate 



Si(.(K)) = -^ 



1 - Re^ ^ — 



2(l + i?e2fi2(K;ap))2 
which is identical to (75) with the appropriate rescaling and rotation. 
In the case of the base flow (A. 13), (A.18) gives the growth rate 
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K(a(K)) = - 



'Re 



1 - 



3Re^ Re^n\K;ap) -1 



4 {l + Re''n^{K;ap)y 



which is identical to (80) after an appropriate rescahng. 

In the case of the base flow (A. 11), since it does not satisfy (A. 14) for c 7^ 1, 
we can not use (A. 18) to find the growth rate. Multiple Scales analysis for 
(A. 11) can be used in a similar way as above and the result is (71). We omit 
this and refer to Manfroi and Young [20] for details. 
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